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Abstract 

The rigorous explanation for the term \t\ 2 @ in the rectilinear diameter equation is given (t = 
(T c — T)/T c , f3 is the critical exponent for the asymptotic form of the equation of state). The 
optimal order parameter, for which the branches of binodal are symmetric is constructed within 
the canonical formalism. It is shown that the ratio of the amplitudes for the diameter singularity 
of the order parameter ~^~ Q! before |t| 1_a and |t| 2 ^, where a determines the behavior of the heat 
capacity and (5 is the critical exponent of the order parameter, takes the universal character modulo 
nonuniversal factor which depends on the thermodynamic class of the corresponding states. The 
analysis of entropy for argon and water leads to j3 = 0.33 and the corresponding amplitude ratio 
» -3.5. 
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I. INTRODUCTION 



The nature of the asymmetry of the coexistence curve in liquids and multi-component 
solutions has a long history, which dates back to the work of Cailletet and Matthias QXJ , where 
the empirical linear law was obtained (see also [0] and references therein). The diameter 
(p^(t) is determined as following: 

^ d \t)= 1 -(^\t) + ^\t)) } (1) 

where yjW(t), i = l,g, are the values of the order parameter on the liquid and gas branches 
of the coexistence curve, considered as a function of the dimensionless temperature t = T ~ Tc 
(T c is the critical temperature). According to the fluctuation theory of critical phenomena 
|| ipW(i) is described by the expansions: 

^ J A* + !%«" + ..., |«|> to _ 

[ D 1 ^ a \t\ 1 ~ a + . . . , \t\<£t G 

where ta is the Ginzburg temperature, and a is the critical exponent for the heat capacity. 
The anomaly, where t = T ~ Tc , of the rectilinear diameter was predicted by Widom 

and Rowlison in Q. Since 1 — a is close to unit (« 0.89), the influence of the singular 
term in (|j) is expected to be essential only in the intermediate vicinity of the critical point. 
However, the careful analysis of experimental data for the rectilinear diameter may display 
another type of the critical singularity: 

V id) (r) = D w \t\ 2 ? + D^ a \r\ l - a + D x r + D 2 r 2 + ..., (3) 

where (3 is the critical exponent for coexistence curve: 2(3 ~ 0.654 < 1 — a (for the review 
see 0). Because of very close values of the exponents of the first three terms one needs 
the data of high precision to get the reliable estimates for the amplitudes Dj. Here we use 
the temperature like variable r, which depends on the reduced temperature t and will be 
specified below. 

The rectilinear diamenter singularity has intensively discussed in H, U in connection with 
the nature of "|t| 2/3 " anomaly and its relation with the Yang- Yang (YY) singularity 0. As 
was shown in || the appearance of the \t\ 2 @ term in the rectilinear diameter demands that 
the equation of state is as following: 

7l (C,r) = n + r 2 - a f ± ^-^^ , r = r(P,T^), C = ({Pi T, /i) , (4) 



where 7r, t and ( are the analytic functions of all three thermodynamic fields P, /i, T. Here 
C is the field conjugated to the order parameter, which in general is not density, r is the 
"temperature "-like variable and 7r is the thermodynamic potential corresponding to them 
(the details see in Thus the scaling law (Q) states that there is no preferred variable 
between P, fi and T ||. Within the approach proposed in |5], || the YY anomaly as well as 
"|t| 2/3 " singularity appear as the result of the nonlinear dependence of the physical quantities 
on the scaling fields due to complete/pressure mixing, which is postulated. However, the 
regular method to construct the variables ( and r was not proposed. 

Here we will not discuss the appearance of the \t\ 2f3 due to nonlinear transformation of 
the laboratory order parameter (of type v = 1/n, where n is the density and v is the specific 
volume). It is clear that e.g. for the Ising model one is free to chose the "non symmetrical" 
order parameter m = m+^m 2 , where m is the magnetization, which gives the corresponding 
r 2/3 t erm for the rectilinear diameter. The choice of the proper order parameter in binary 



mixtures which restores the symmetry of the binodal was discussed in |T0| , |TT|| , though any 
procedure for constructing such an order parameter was not given. The problem is to find 
such a microscopic observable r\ so that: 

along the coexistence curve. Then in the framework of the RG theory the asymptotic 
equation of state is fully symmetrical: 

(,) = ±r" ff ,(A_) (5) 

where h and r are the external field and temperature like scaling fields correspondingly. 
However the direct comparison of the theoretical prediction with the experimental data is 
impossible since the explicit form of the order parameter r] and the fields h and r remain to 
be unstated. 

The symmetrical form of the Eq. (H) is based on the symmetry property of the effective 
Hamiltoniian H,[(p(r)] = 7~C[— y?(i")] along the binodal. In particular the Landau- Ginzburg 
Hamiltonian satisfies this condition ||. 

In this paper we give the rigorous formalism for the construction of the scaling fields r and 
h. Within such a formalism the nonlinear mixing of the laboratory fields arises naturally. 
In particular we show that both \t\ 2f3 and \t\ l ~ a anomalies are of the same nature. 



Further we demonstrate this statement using the entropy as the order parameter. It is 
natural that the entropy is more sensitive to the asymmetry in the "particle- hole" configu- 
ration than the density. The latter is just the average of 1-particle distribution function for 
the particles. 

The canonical transformation ip — > r](ip) of the initial (laboratory) order parameter <p is 
the central point of the proposed approach. Due to this the local part of the fluctuation 
Hamiltonian of liquids reduces to the canonical form, which is identical to the local part 
of the Landau- Ginzburg (LG) Hamiltonian. In canonical variables the coexistence curve is 
fully symmetrical, as it takes place for the Ising model. The non-zero values of the rectilin- 
ear diameter arise as the result of the inverse transformation 77 — > ip to the laboratory order 
parameter. In the framework of the canonical formalism the isomorphism of critical fluctu- 
ations in liquids and the Ising model is manifested in the similarity of their thermodynamic 
potentials as well as their Hamiltonians. Within the perturbative approach it can be shown 
that the asymmetry in main order depends on the value of the coefficient 05 at ip 5 term in 



the initial Hamiltonian. This is accordance with the results of [12 



The construction of the canonical order parameter r/ for which rf^ = is discussed below. 



II. CANONICAL FORMALISM 



The Hamiltonian for the real fluids can be represented as following: 

H[p(r)} = J ( hMv)) + h ql (p(r)) ) dV . (6) 

where 

00 

n 

n=l 

is the local density and 

00 

' n 

n=l 

is the quasilocal part in small gradient approximation. The Hamiltonian ([D includes odd 
power terms of the order parameter responsible for the "asymmetry effect" in the critical 
behavior of the system. It should be noted that there are two kinds of asymmetry terms in 
the Hamiltonian. They are the local and the quasilocal terms of odd powers. The importance 
of the lasts for the RDS was shown in |12|, |13|| . 
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Note that the local term ip 5 leads to the mixing of the initial thermodynamic variables [l 
and T. The quasilocal term <£>(Vy?) 2 leads to the I — a singularity of the rectilinear diameter, 
though 2(3 anomaly was not obtained there. 

Now let us appeal to the ideology of the Catastrophe Theory. It gives the ground to 
reduce the initial Hamiltonian given by the infinite series Eq. (Q) to the canonical form. 
For this purpose the order parameter should be subjected to the so called local canonical 
transformation 

V (v) = <p(r) + I T 2 ^ 2 (r) + i T 3 ^ 3 (r) .... (7) 
As a result the local density of the Hamiltonian takes the form: 

h!r\rj) = A lV (r) + ^rf{v) + ^ V \r) . (8) 

with A 4 = const. Using the procedure described in [|TJJ one can show that 

r 2 =A X - a x , 

r 3 = ^ 2 -iA 2 + ^a 1 r 2 -ia 2 + ia 2 + ^r 2 , (9) 



1 15 
T 4 = — Ax A 2 + - A 3 , - - ax T 2 

2 6 1 2 2 



4 11 
+ a 2 T 2 - a 2 r 2 + - ax T 3 - - a 3 + - al a 2 

-hal- -r| + -r 3 r 2 , (10) 

6 1 2 2 3 ' v 1 

In this paper we illustrate the explicit construction of the transformation (0) which reduce 
only the local part of the Hamiltonian to the canonical form. In such a case the coefficients 
Ax and A 2 should vanish simultaneously with ax and a 2 : 

Ax = 0, A 2 = 0^ax = 0, a 2 = . (11) 

These conditions correspond to the invariance of the locus of the CP in mean-field approxi- 
mation. 

To find the explicit expressions for the coefficients Ax, A 2 and A 4 we demand the fulfill- 
ment of the following conditions: a) the conservation of the local partition function: 

00 00 

J dxexp (-/i£ n) (z)) = J dyexp(-h loc (y)) . (12) 



b) one-to-one correspondence (bijectivity) of the canonical transformation (0): 

77 <p 



J exp (-h[2 n \z)J dz = J exp (-h t (z)) dz . 



(13) 

0, which in the 



c) the condition of the invariance of the phase coexistence line f(a±, 02) 
canonical variables is represented by the equation A\ = 0. 

Note that similar idea of transformation of the variable reducing the density distribu- 
tion function to simpler (gaussian) form was used in |I5|. In fact the existence of such a 
transformation is guaranteed by the Radon-Nikodim theorem l£j . 



Basing on these conditions one can show (see some details in JTj]) that the coefficients of 
the canonical form satisfy the relations: as: 



A, 



7T 



r(3/4) / dyexp - £ 



n=4 



4 ' 



(14) 



IA2 A T . 
— e 2A i Ki 



2A d 



exp (— hi(z)) dz 



(15) 



/(ai,a 2 )=0 



Here K n is the McDonald function of ra-th kind and T(x) is the Euler gamma function |17 



Neglecting the variation of nonlocal part for the fluctuation part of the Hamiltonian and 
after the rescaling we get the Landau- Ginzburg (LG) form: 



n LG W) 



A xV {v) + ^-r?(v) + ^ V \v) + i (V^(r)) 2 ) dV 



(16) 



equivalent to the Hamiltonian of the Ising model. Thus the possibility to transform the 
Hamiltonian to the canonical form (§) is equivalent to the isomorphism of the critical be- 
havior of liquids with that for the Ising model. Usually, for the liquids the Hamiltonian @ 
is obtained via grand canonical ensemble, which means that the coefficients a, and Ai are 



analytic functions of the thermodynamic variables fi and T fl8 |. 

Applying the standard scaling arguments (see || |19|) to the canonical Hamiltonian 
we get the fluctuation contribution to the thermodynamic potential: 

A 1 



$(A 1 ,A 2 )=$ + \A 2 \'- a f, 



\Ao\P+"> 



(17) 



Here $ is the thermodynamic potential which corresponds to the canonical coordinates 
A±, A2, which can be considered as the generalized external field and the temperature con- 
jugated to the canonical order parameter r\. The potential $ does not coincide with the 



pressure since the order parameter 77 conjugated to it is not the density. The potential $ is 
rather the function of all three variables P,T,fi and corresponds to Eq. (^). Of course, one 
can choose another set of variables, e.g. P, T. It is in accordance with the P, /1, T invariance 
hypothesis mentioned above and lead as was shown in M to the "t 2 ^ " anomaly in the recti- 
linear diameter of the density and as well as the entropy. We show below that the canonical 
formalism explicitly shows that both 2/3 and 1 — a singularities of the rectilinear diameter 
of the density as the order parameter are of the same nature. They are generated by the 
asymmetry of the initial Hamiltonian. 

Note that the local transformation generates also the quasilocal term r](Vi]) 2 multiplied 
by the coefficient which is proportional to T 2 . In this case all local odd power terms as 
well as the term ^(Vr/) 2 can be canceled. The coefficient A 4 is determined by the last 
condition. It means that the condition of the invariance of the mean-field critical point 
Eq. ([11]) is not assumed. The coefficients Ai and A 2 are correspondingly modified, which 
lead to the shift of the mean-field critical point. The appearance of such shift is connected 
with the inclusion of the quasilocal interaction between modes of the order parameter. It 
is necessary to emphasize that this shift arises before the renormalization procedure. The 
detailed analysis will be the subject of separate work. 



III. CANONICAL FORM FOR VAN DER WAALS EOS 

Let us obtain the explicit expressions for the coefficients of the canonical form of the 
Hamiltonian for a case, when 1) the density of the local Hamiltonian has the truncated 
9? 6 -form: 



n 

n=l 



2) the coefficients a n are taken from the van der Waals equation of state: 

= dl73- 8^ (18) 
where p, v and t are the dimensionless pressure, the specific volume and the temperature 

p = P/P c , t = T/T c , v = V/V c 

reduced to the coordinate of the vdW critical point (see e.g. [p0|). It can be shown that: 

_p-l + 4r r t _31 + 9r 

a i — — T~ j a 2 — -J 7- - — , a 3 — -J7- — , a 4 — -— 

1+r 1+r 1+r 21+r 
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21 1 + §r _ 99 1 + §r 

Applying the procedure described above we get: 



(o) 

A A « af> + 1.35 —f= » 15.2 (19) 



(0) 



where is the value of the coefficient at the critical point. Other coefficients are as 



following: 



2 a 3 2 4 0) „ , , 

A x « fll + — ^= + - , A 2 « a 2 . 20 

3 v^ . /lo) 5 a (°) 



From Eq. © and Eq. © it follows that 



r 2 «-aS 0) 1 + I (21) 

5 5 3v^F f 0)i ri 1 



a 5 V a 4 



It is necessary to emphasize that these estimations are valid only in the mean field approx- 
imation. For the accurate account of the fluctuational effects the quasilocal part of the 
Hamiltonian should be also taken into consideration Wl . 



IV. THE 2/3 RECTILINEAR DIAMETER SINGULARITY WITHIN THE 
CANONICAL FORMALISM 

The coexistence curve determined by Ai([i, T) = in the variables (Ai, A 2 , rj) is naturally 
symmetrical. The canonical order parameter average value ( r\ ) in accordance with Eq. fllj 
and the thermodynamics is: 
<9$ 



(v) = - 



dA t 



±177*0.(0)+..., f=A 2 \ Ai=0 = at + o(t), (22) 



.41=0 



where g s {%) = f' a {%)- The function f s (x) possesses obvious symmetry f s (x) = f s (—x), which 
means that in the canonical variables the "particle-hole " symmetry is restored. The situation 
changes if we return to the initial order parameter: 

1 



<P = V~ 2 r 2 ^ + 



Since r\ is the symmetrical order parameter, for the rectilinear diameter ^ > we get: 

^ d) = -\^W)+0{{^)) . (23) 



From Eq. (|23|) and the short distance behavior of the correlator G(r) = (r](r) rj(0)) — 
(r](r)) (77(0)) for the order parameter of the LG Hamiltonian (||) which determines the 
critical behavior of the fluctuation {rf ) — ( rf) 2 (see e.g. 0, j21f): 

< V 2 > - ( V) 2 = G(r -> 0) = J G(k) dk ~ Z s (0) A\ a + less singular terms (24) 

where l s (0) is some constant. Therefore, along the coexistence curve we obtain: 

^(d) = _lf 2 [^(o)|f| 2/3 + Z fl (0)|f I 1 ""] + less singular terms, (25) 

where r 2 the value of corresponding quantity along the coexistence curve: 

f 2 = r 2 U l=0 . (26) 

Thus the amplitudes Di are: 

= -if 8 Z.(0), ^ = -if 2 ^ 2 (0). (27) 

and either both zero or nonzeroth. Note that within such an approach both f 1_a and f 2/3 
anomalies contribute to the rectilinear diameter for any parameter of state tp chosen as the 
initial order parameter, e.g. entropy, density, dielectric permittivity etc. Thus we state that 
these anomalous terms are of the same general nature and should be treated simultaneously 
for the "nondegenerate " systems in terms of the Catastrophe Theory. They appear as the 
defect of the improper choice of the order parameter and are generated by the asymmetrical 
part of the Hamiltonian. This is different from the result of work || where t 2 ^ anomaly is 
due to nonlinear mixing while anomaly is due to the linear mixing of the fields. 

It is important to note that all results discussed in this Section do not depend on the 
method for the determination of the coefficient IV In general case T 2 determined by the 
contributions generated by both local and quasilocal terms in the initial Hamiltonian. 

Note that since l s (0) < (see 0), the amplitudes Di_ a and D 2 p are of opposite signs. It 
is important that their ratio is universal: 

^ = M^<0. (28) 
D20 9, 2 (0) 

Processing the experimental data is usually performed in terms of the reduced temperature 
variable t, i.e. in Eq. (p5|) the change f — > t is made. Therefore the ratio of corresponding 
amplitudes 

Dt a = A-a a x ~\ D^ = D 2(3 a 2 ^ (29) 



is non universal because of Eq. ( pq) . Therefore the ratio is as following: 

)(*) 

(30) 



D *- a „ 1.(0) 



D% 9 2 s(0) ' 

where the nonuniversal factor a is given by Eq. (|22|). The values of the amplitudes D 2 /3 and 
Di- a depend on T2- The latter can be either positive or negative depending on the details of 
the intermolecular interactions. This may explain the fact that for molecular liquids D 2 p > 
while in the case of liquid metals the opposite sign takes place ]22|, Q . This may attributed 



to the difference in the short range interaction in these systems. For the molecular fluids 
the pure hard core interaction takes place while in liquid metals the polarizability effects are 
essential due to incomplete electronic shells e.g. in liquid Hg p4|, p5fl . Strong polarizational 
effects in liquid metals and ionic fluids give rise to screening and the formation of bound 
states (dimers, trimers etc.). In such a situation one can expect that the quasilocal part of 
the Hamiltonian for such systems differs very much from that of molecular liquids because 
of strong dispersive interparticle interactions. 

Though the functional form f(t) is non universal, it is the same for the systems, which 
obey the law of corresponding states. This inference is in correspondence with the estimates 



made in for molecular and ionic liquids for which the density as the order parameter 
was used: 

%^«-10. (31) 

D2/3 

V. ANALYSIS OF THE DATA FOR ARGON AND WATER 

To confirm the results obtained we have analyzed the experimental data on the entropy 
s for water and its saturated vapor (see Fig. 0). In this case the initial order parameter 
ip is naturally defined as <p = It is natural that the entropy is more sensitive to 

the configuration "hole-particle" than the density, which is just the average of 1-particle 
distribution function. 

The evaporation heat q and rectilinear diameter are determined in the standard way: 
q = T c s c (tf^ 9 ' — ip( l >) and Sd(t) = ^(S^ + S®). The behavior of ip^ for water and argon 
is presented in Fig. |I[ In accordance with said above the corresponding experimental data 
were fitted by the formulas: 

\t\P + q 2 \t\P +A + q 3 t + o(t) (32) 



1 

t7 = (/i 
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. 1 




. 1 



FIG. 1: The rectilinear diameter of entropy for water (solid) and argon (points) according to data 



of H 



and 



S-W(t) = sg|t| 2/3 + Sfljtl 1 -" + S[ t] t + o(t) . 



(33) 



Note, that as follows from the representation: 



S = c v In T + f{n) 



(34) 



Si_ a < 0. In contrast to the coefficients -Di_ a which is generated by the asymmetry of 
the Hamiltonian, the coefficient Si_ a does not vanish in symmetrical case, e.g. in the Ising 
model, and correspondingly: 



c(*) _ q{ s V m ) I q(asym) 

J l-a ~ °l-a ~r °l-a 



(35) 



It is obvious that the symmetrical part S^™ is generated by the heat capacity singularity. 
The asymmetrical part of the amplitude is of the same structure as Eq. (EW) with corre- 
sponding ratio (j30|). At that we put a and A to be known: a = 0.11 and A = 0.5. The 
analysis shows that (3 is the essential function of the temperature interval t exp , inside which 
the experimental data were taken into account. The fitting was performed for (3 and the 
coefficients <Sj in order to achieve the minimum for the standard deviation. The value 
of such a minimum depends on the interval of interpolation. Both for argon and water the 
optimal interval corresponded to < t < £* with t* ~ 0.2 with close values of coefficients 
though the data set for argon is smaller than that of water. The ratio (|30"D both for argon 
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and water takes close values: 

— ^ & -3.4 , -^«-3.7, (36) 

which is in correspondence with the result (|28|). The difference is attributed to the nonuni- 
versal factor a (see Eq. (^6])) which depends on the details of the canonical transformation. 
In order to distinguish between the t 2 * and t l ~ a terms for rectilinear diameter of the den- 
sity we subtracted the analytic background terms taking into account the data beyond the 
region t > to- Such a subtracting allow to fix approximately the value to because the linear 
analytic term is important even in fluctuational region. We fit the data in the region t > t* 
with the Taylor expansion in t up to t 5 order and check the difference between this regular 
expansion and the data. At t* > 0.2 no significant deviation between them was detected. 
Therefore we can conclude that the crossover takes place at to ~ 0.2. 

It should be noted that it is impossible to fit the data with the the singular t 1_a! -term 
only. The t 2/3 -term is essential and together with the analytical terms allows to fit the data 
within rather broad interval < t < 0.2 correctly both for the argon and water. Note that 



in approach used |27|> £3] in the Ginzburg temperature depends on the choice of the order 
parameter and may differ substantially depending on the choice of the truncated expansion 
for the initial Hamiltonian. 

Note that the existence of additional t 2/3 term is also implicitly assumed in empirical 
equation of state for water of Wagner&Pruss, which fits the data in whole coexistence region 
up to the triple point with the value j3 = 0.33, though the term ~ was not taken into 
account ||26|| . From this point of view the equations of state in broad interval, which uses 
the crossover functions are preferable |p7| , [29| . 

The authors cordially thank Prof. M. A. Anisimov and Prof. H. E. Stanley for the 
discussion of obtained results. 
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